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ABSTRACT 

We present a redshift z = 6.535 galaxy discovered by its Lyman-a emission in 
a 9180A narrowband image from the Large Area Lyman Alpha (LALA) survey. The 
Lyman-a line luminosity (1.1 x 10 43 erg s _1 ) is among the largest known for star forming 
galaxies at z ~ 6.5. The line shows the distinct asymmetry that is characteristic of 
high-redshift Lyman-a. The 2a lower bound on the observer-frame equivalent width 
is > 530 A. This is hard to reconcile with a neutral intergalactic medium unless the 
Lyman-a line is intrinsically strong and is emitted from its host galaxy with an intrinsic 
Doppler shift of several hundred km s _1 . If the IGM is ionized, it corresponds to a rest 
frame equivalent width > 40A after correcting for Lyman-a forest absorption. We 
also present complete spectroscopic followup of the remaining candidates with line flux 
> 2 x 10~ 17 erg cm -2 s~ 1 in our 1200D' narrowband image. These include another galaxy 
with a strong emission line at 9136A and no detected continuum flux, which however is 
most likely an [O III] A5007 source at z = 0.824 based on a weak detection of the [O III] 
A4959 line. 

Subject headings: galaxies: high-redshift — galaxies: formation — galaxies: evolution 
- cosmology: observations — early universe 
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1. Introduction 

Observational study of the redshift range 6 < z < 7 is crucial for understanding the reionization 
of intergalactic hydrogen, which was the most recent major phase transition for most of the baryonic 
matter in the universe. While polarization of the microwave background measured by the Wilkinson 
Microwave Anisotropy Probe (WMAP) satellite indicates that substantial ionization had begun as 
early as z ~ 15 (Spergel et al. 2003), both the opaque Gunn- Peterson troughs observed in the 
spectra of z > 6.3 quasars (Becker et al. 2001; Fan et al. 2002) and the high temperature of the 
intergalactic medium at z ~ 4 (Theuns et al. 2002; Hui & Haiman 2003) imply that a large part 
of the reionization was relatively recent, with substantial neutral gas lasting up to z ~ 6. 

Lyman-a emission has proven the most effective tool so far for identifying star forming galaxies 
at redshifts z > 6; indeed, all but two of the ~ 6 galaxies that have been spectroscopically confirmed 
at z > 6 were found through their Lyman-a line emission in either narrowband or spectroscopic 
searches (Hu et al 2002 [H02]; Kodaira et al 2003 [K03]; Cuby et al 2003; this work; Kneib et al 
2004), and at least one of the continuum-selected sources is also a Lyman-a emitter (Cuby et al 
2003). 

We present the extension of the Large Area Lyman Alpha (LALA) survey to the z « 6.5 
window. We have obtained spectra for each of the three z ~ 6.5 candidates that pass all photometric 
selection criteria. Two of these show strong emission lines in the narrowband filter bandpass. In 
one case, the line is identified as Lyman-a at z ~ 6.535 based on its asymmetric profile and the 
absence of other detected lines down to faint flux levels. The second object is identified as an [O III] 
A5007 source at z ~ 0.824 based on a probable (4a) detection of the [O III] A4959 line. 

2. Imaging Observations and Analysis 

The LALA survey s z ~ 6.5 search is based on a deep 9180A narrowband image of our Bootes 
field, which is 36' x 36' with center at 14:25:57 +35:32 (J2000.0). The image was obtained using the 
CCD Mosaic- 1 Camera at the 4m Mayall Telescope of the Kitt Peak National Observatory together 
with a custom narrowband filter, whose central wavelength of A c ~ 9182A and full width at half 
maximum (FWHM) of 84A place it in a gap between the night sky airglow lines. 

Narrowband imaging data were obtained on UT 2001 June 15-17 and 24-25 and 2002 April 5 
and 18-19. The total exposure time was 28 hours. 
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Data reduction followed the same procedures used in the z = 4.5 and z = 5.7 LALA images 
(Rhoads et al. 2000; Rhoads & Malhotra 2001). To summarize, we remove electronic crosstalk 
between Mosaic chip pairs sharing readout electronics; subtract overscan and bias frame corrections; 
and flatfield with dome flats. Next, we subtract a pupil ghost image caused by internal reflections 
in the KPNO 4m corrector. We remove residual, large scale imperfections in the sky flatness using 
a smoothed supersky flat derived from the science data, followed by subtraction of a polynomial 
surface fit to the sky flux. A fringe frame was constructed from the data, fitted to the blank sky 
regions of each exposure, and subtracted. Cosmic rays were rejected in each exposure using the 
algorithm of Rhoads (2000). 

The world coordinate systems of individual frames were adjusted to the USNO-A2 star catalog 
(Monet et al 1998). Satellite trails were flagged by hand for exclusion from the final stacks. The 
eight-chip Mosaic-1 images were then mapped onto a rectified, common coordinate grid using the 
"mscimage" task from the IRAF MSCRED package (Valdes & Tody 1998; Valdes 1998). The 
seeing, throughput, and sky brightness in each exposure were measured and used to compute a set 
of weights for image stacking using the "ATTWEIGHT" algorithm (Fischer & Kochanski 1994). 
Finally, the exposures were stacked, with an additional sigma clipping applied to reject discrepant 
data not previously flagged. The majority of the weight comes from the final two nights' imaging 
data, and only 24% of the total weight is from the 2001 observing season. 

The final stack has seeing 1.02" and a 5<r sensitivity limit of 0.7//Jy (measured within a 9 pixel 
[2.32"] diameter aperture), corresponding to 2.0 x 10~ 17 erg cm -2 s _1 for a pure line source at these 
wavelengths or to an AB magnitude of 24.3 for a pure continuum source. The magnitude zero point 
was determined by comparison with a deep z' filter image obtained earlier (see Rhoads & Malhotra 
2001), which was in turn calibrated to z' filter standard stars from SDSS standard star lists. 

3. Candidate Selection 

Lyman-a galaxy candidates were selected following the criteria used successfully in lower red- 
shift LALA searches (Rhoads &; Malhotra 2001; Dawson et al. 2004). In summary, these include: 
1) Significance of the narrowband detection > 5a; 2) A narrowband excess of at least 0.75 mag- 
nitude, so that > 50% of the narrowband flux comes from an emission line; 3) Significance of the 
narrowband excess > 4a; and 4) No detection in filters blueward of the expected Lyman break 
location at the > 2a level. 

To implement the blue flux "veto" criterion, we used a weighted sum of six NOAO Deep Wide- 
Field Survey (NDWFS; Jannuzi k, Dey 1999) and LALA survey filters blueward of 6850A, where 
the Lyman break would fall for galaxies whose Lyman-a line falls within the 9180A narrowband. 
The weights were chosen to achieve optimal depth for objects of approximately median color. Most 
of the weight in this combined "veto image" comes from the NDWFS Bjy filter, with additional 
contributions from NDWFS R band data and the LALA survey's V band, and three of the redshifted 



-4 - 



Ha narrowband filter observations. This stack reaches a final 2a limiting AB magnitude of 27.1 at 
an effective central wavelength near 5000A. Formally, some flux may be expected in the R band 
filter for z ~ 6.5 sources, but in practice, intergalactic hydrogen absorption attenuates this flux so 
severely that it is quite safe to include R band data in the veto stack as we have done. 

Photometry for candidate selection was done using SExtractor (Bertin & Arnouts 1996). All 
color tests were applied using a 9 pixel diameter aperture, corresponding to 2.32", and using 
SExtractor's two-image mode to ensure photometry from the same regions in all filters. 

In addition to high redshift Lyman-a sources, three types of contaminants may enter the 
sample. The first is noise spikes in the narrowband image. Assuming Gaussian statistics, we 
should expect ~ 1 noise peak above our 5a detection threshold among the ~ 10 65 independent 
resolution elements in the narrowband image, and more may be found if the noise properties of 
the image are not perfectly Gaussian. The second is time variable sources (either transient or 
moving objects). These may mimic emission line objects in our catalog because different filters 
were observed at widely varying times. For typical LALA survey depths, the expected rate of high 
redshift supernovae will be ~ 1 per filter per field (Riess, personal communication; see also Strolger 
et al 2004). The third is low-redshift emission line galaxies of extremely high equivalent width. 
The expected numbers of such foreground objects are difficult to estimate from either theory or 
present data — the LALA data themselves are likely to considerably refine the statistics of such 
contaminating sources. 

These criteria yielded a total of three high quality candidate z ~ 6.5 Lyman-a galaxies. For- 
mally, six objects passed selection using an old veto filter stack from 2001, but two of these were 
ruled out completely and a third strongly disfavored by detections in a new NDWFS R band image 
from 10 April 2002. One of these excluded candidates is almost certainly a supernova at moderate 
redshift; such an event would not be surprising given the area and depth of our survey. Properties 
of the four best candidates are summarized in table 1. 
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Table 1. Photometry 



ID 


918nm flux 


z' flux 


/ flux 


R flux 


Photometric 


W r A cst (2cr) 


Nature 




(A*Jy) 






(^Jy) 


line flux a 


limit (A) 




LALA J142441.20 


0.763 


-0.214 


-0.0145 


-0.0372 


4 to 6.6 b 


> 620 


2 = 0.824 


+353405.1 


±0.139 


±0.111 


±0.045 


±0.029 








LALA J142442.24 


0.767 


-0.081 


-0.064 


-0.061 


2.26 


> 70 


z = 6.535 


+353400.2 


±0.137 


±0.113 


±0.045 


±0.028 


±0.40 






LALA J142544.41 


0.706 


-0.051 


-0.073 


-0.021 


2.08 




Noise spike 


+353327.7 


±0.138 


±0.106 


±0.044 


±0.029 


±0.41 






LALA J142610.55 


0.797 


-0.137 


0.011 


0.057 


2.35 




Unknown 


+354557.6 


±0.121 


±0.110 


±0.046 


±0.030 


±0.36 







a Units of 10~ 17 erg cm -2 s" 1 . 

b The line for LALA J142441.20+353405.1 falls on the edge of the narrowband filter, where throughput 
is a rapid function of wavelength. The dominant uncertainty in the photometric line flux is the resulting 
throughput correction at 9136A, and we quote a range of flux values based on the plausible range of this 
throughput correction. See text. 

Note. — Photometric properties of the candidate z w 6.5 LALA sources. The 2a lower bounds on equiv- 
alent width are derived from the narrow- and broad-band photometry for LALA J142442. 24+353400. 2, and 
from DEIMOS spectra for LALA J142441. 20+353405.1, and are corrected to the spectroscopically determined 
redshift. Tabulated values are not corrected for IGM absorption. 
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4. Keck and Gemini Spectra 

Two of the z=6.5 candidates (LALA J142441. 20+353405.1 and LALA J142544.41+353327.7) 
were observed using the DEIMOS spectrograph (Faber et al. 2002) on the Keck II telescope of the 
W. M. Keck Observatory on the night of U.T. 2003 May 01. Observations of LALA J142441. 20+353405.1 
were also obtained on U.T. 2003 Apr 01. All observations were made using the 600 1/mm (A& = 
8500A) grating through slitmasks with slit widths of l.Oarcsec. This setup gives a dispersion of 
~ 0.63A/pixel and a spatial plate scale of 0.1185arcsec/pixel. The total exposure times were 
2.5 hours for LALA J142441. 20+353405.1 (1 hour on Apr 01 and 1.5 on May 01) and 2.0 hours 
for LALA J142544.41+353327.7. These were broken into individual exposures of 1800s, with 2.5 
arcsec spatial offsets between exposures to facilitate the removal of fringing at long wavelengths. 
The nights were clear and the seeing was typically 0.5-0.8arcsec. The data were reduced using the 
U.C. Berkeley pipeline reduction software (Davis et al. 2004) adapted from programs developed 
for the SDSS (Buries and Schlegel 2004), and further adapted to our observing mode. The spectra 
were extracted and analyzed using IRAF (Tody 1993). 

The candidates LALA J142441. 20+353405.1 and LALA J142442. 24+353400. 2 were observed 
on U.T. 2003 May 8 and June 29 using the GMOS spectrograph on Gemini North (Hook et al 
2003) in Nod-and-Shufhe mode (Glazebrook & Bland-Hawthorn 2001). The observations were 
made using the R400+G5305 grating (400 lines/mm, central wavelengths at 756 - 768 nm) through 
slitmasks with slit widths of l.Oarcsec. On each night, the total exposure time was broken into 5 
individual exposures of 900s per nod position. Subsequent exposures were offset by 0.2 arcsec in 
the spatial direction and 3 nm in wavelength, in order to remove the instrumental features (such 
as the horizontal stripes in the individual exposures) and to cover the chip gaps. Since the two 
objects appear on both the target and sky exposures (i.e., on both shuffle A and shuffle B), the 
effective exposure time for the two sources are 9000s or 2.5 hours per night, for a grand total of 5 
hours' on-source integration. 

Additionally, a second GMOS slit mask was observed for a total of 3.5 hours total integration 
time on UT 2003 July 02 (0.5 hour) and July 03 (3.0 hours), covering the (lower grade) candidate 
LALA J142610.55+354557.6. 

The data were reduced using the Gemini GMOS packages (vl.4) within IRAF through the 
standard procedures (http://www.gemini.edu/). The wavelength calibrations were performed based 
on the CuAr lamp spectra and were double checked against the night sky lines. The uncertainties 
in the extracted spectra were obtained empirically by measuring the RMS flux among all spatially 
distinct pixels at each wavelength in the rectified 2-D spectra. 

Additional targets were of course included on all slit masks, including lower redshift emission 
line galaxies, continuum-selected Lyman break galaxies and intermediate redshift elliptical galaxies, 
and X-ray sources from the 172 ksec LALA Bootes field Chandra ACIS observation (Wang et al. 
2003; Malhotra et al. 2003). Results for these other targets will be presented elsewhere. 
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5. Spectroscopic Results 

Two of our candidates, LALA J142442. 24+353400. 2 and LALA J142441. 20+353405.1, yielded 
emission lines within the bandpass of the narrowband filter. Of these, one 
(LALA J142442. 24+353400. 2) is confirmed as a Lyman-a galaxy, while the other 
(LALA J142441.20+353405.1) is identified as [O III] A5007. 

LALA J142442. 24+353400. 2: LALA J142442. 24+353400. 2 was confirmed by Gemini-N + 
GMOS on 2003 May 8 and 2003 June 29. No DEIMOS data were obtained for this object. There 
is a single isolated emission line at 9160A, corresponding to z = 6.535 (see figures 1, 2). No other 
lines are seen in the observed wavelength range, 6098A < A < 10484A. 

There is no evidence for continuum emission in either our images or GMOS spectrum. The 
observer-frame equivalent width of the line is > 530A at the 2a level. We obtained this limit our 
LALA survey narrow- and broad-band photometry as follows. First, let and t\ tTl be the system 
throughput in the broad and narrowband filters as a function of wavelength, which we base on 
the filter transmission curves and quantum efficiency of the Mosaic camera. Then given the line 
wavelength A^, we define Wf, = (/ t\^d\) /t\^(\e) and w n = (J t\ jn dX)/t\ jn (\£). For a narrow line 
atop a flat (i.e., constant fx) continuum, one can show that EW = w\>w n (f\^ n — f\,b)/(f\,b w b ~ 
f\,nW n ), where and f\ >n are the flux densities measured in the broad and narrow band filters. 
We compute a 2a limit on EW by inserting f\j, — > f\ :Zj0 bs + 2df\ >z into the above expression. We 
ignore the uncertainty in the narrowband flux f\ :Tl , because it is too small to much affect the EW 
measurement. For Lyman-a at z = 6.535, the rest frame equivalent width would be > 70A(2cr) 
if we do not correct for attenuation by the intergalactic medium (IGM). The IGM is expected to 
attenuate the z' filter flux by a factor of ~ 0.36 at z = 6.535 (Madau 1995). The correction to the 
line flux depends on the IGM neutral content, IGM velocity structure, and intrinsic line profile. 
For a symmetric line centered on a galaxy's systemic velocity, the line flux would be reduced by a 
factor of 0.51 in a predominantly ionized IGM (i.e. all of the red side flux but only ~ 2% of the blue 
side flux would be transmitted). The IGM-corrected rest equivalent width would then be > 40A 
(2a), which we regard as a conservative lower bound. A largely neutral IGM or a line profile not 
symmetric about the systemic velocity would change this result appreciably (see section 6.2). 

The line width in this source could in principle be consistent with [O II] AA3726,3729. However, 
this would imply a rest frame equivalent width > 2 15 A (2a), which is outside the usual range for 
this line: two large continuum-selected samples (Cowie et al. 1996, Hogg et al. 1998) and one Ha 
line selected sample (Gallego et al. 1996) found no [O II] AA3726,3729 sources with EW > 140A 
and a very small minority with 100 < EW < 140A. Additionally, star forming galaxies usually 
show blue continuum emission, so that our and R filter nondetections again argue against [O II] 
AA3726,3729 models. 
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Fig. 1— Spectrum of LALA J142442. 24+353400. 2, obtained with the GMOS spectrograph on 
Gemini North. The solid histogram shows the measured flux; the dotted line shows the one-sigma 
flux uncertainties. The asymmetry of the line is clearly visible. There is no emission at 9072A, 
where the [O III] A4959 line would be expected if the 9160A feature were [O III] A5007. The 
conversion from counts to physical units was determined by requiring the integrated line flux to 
match that determined from our narrowband photometry. 
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Fig. 2.— 2 dimensional GMOS spectrum of LALA J142442. 24+353400.2. Night sky lines are well 
subtracted here, but their locations are plainly visible as bands of enhanced photon noise. Arrows 
mark the locations where the H/3 and [O III] A4959 lines would be expected if the primary line at 
9160A were [O III] A5007. 



-10- 



The LALA J142442. 24+353400.2 emission line has a measured width of 13A FWHM. This 
substantially exceeds the instrumental resolution (estimated at 6A based on the GMOS spectrum 
of LALA J142441. 20+353405.1), and thus allows quantitative asymmetry measurements. We apply 
the line asymmetry statistics presented in Rhoads et al (2003). First among these is C sym , the 
likelihood of obtaining the observed data under the best fitting symmetric line model. We find 
C sym « 0.02 for LALA J142442. 24+353400. 2, supporting the visual impression of asymmetry. 

To go further and quantify this asymmetry, we have further developed the asymmetry statistics 
a\ and a/ that were presented in Rhoads et al. (2003). These statistics depend on determination of 
the wavelengths X p where the line flux peaks and Aio,&, Aio, r where the flux drops to 10% of its peak 
value. To mitigate the effects of noise in the spectrum, all three wavelengths can be measured after 
a light smoothing. The asymmetry statistics are then defined as a\ = (Aio, r — A P )/(A P — Aio^) and 
af = fxp°' r f\d\/ f£* o b f\d\. To determine the appropriate smoothing, we simulated observations 
of a line having a truncated gaussian profile (described below) and the (wavelength-dependent) 
noise level achieved in our GMOS data. We find that the asymmetry of the line is most reliably 
recovered using a Gaussian smoothing near 3 pixels (FWHM). We also experimented with the 
threshold level for the Aft and A r measurements, but found no clear benefit to raising this threshold. 
Using 3 pixel smoothing, we find a\ = 1.78 ± 0.71 and af = 1.54 ± 0.53, where the error bars are 
drawn from simulated spectra and defined as a "percentile-based sigma," i.e., 84% of simulated 
spectra lie within ±5a of the median simulated a statistic. 

The LALA J142442. 24+353400. 2 line can be well fit with a truncated gaussian model of the 
form suggested by Hu et al. (2004): The profile is taken to be a Gaussian of width a\ on the 
red side of a peak wavelength Ao and a step function to zero intensity blueward of the peak, 
but then convolved with a second Gaussian of width 02 to reflect instrumental resolution ef- 
fects. In LALA J142442. 24+353400. 2, we obtain x 2 /(degree of freedom) = 0.93 using o x = 8k, 
o"2 = 1.2A, and Ao = 9158A. If we fix 02 = 2.5A (the value obtained by assuming the emission 
line in LALA J142441. 20+353405.1 is unresolved by GMOS), the best fit deteriorates slightly to 
X 2 /(degree of freedom) = 1.17. 

Given that we see no other significant lines in the spectrum of LALA J142442. 24+353400. 2, 
we can rule out [O III] A5007 and Ha interpretations quite strongly. The remaining possibility, 
besides Lyman-a, is [O II] AA3726,3729. However, the line ratio for [O II] AA3726,3729 sources 
under usual astrophysical conditions is fows/ fz72& ~ 1-3, yielding an asymmetry opposite that of 
Lyman-a when the doublet is blended. We have tried modelling the data with [O II] AA3726,3729 
doublets, broadened by the instrumental resolution and an internal velocity dispersion (taken as a 
free parameter). Fixing the line ratio at 1.3 makes it impossible to achieve % 2 /d.o.f. < 1.9, which 
is clearly worse than the truncated Gaussian model for Lyman-a. If we allow the line ratio to vary 
arbitrarily, we can achieve x 2 /d.o.f. = 1.0, but only for 73729/73726 ~ 0.4, which would require an 
electron density > 10 4 cm~ 3 * (T/le4)~ 1 / 5 (e.g., Keenan et al 1999). Extragalactic HII regions do 
not approach these densities, and their line ratios cluster in a tight range 1.2 < Z3729//3726 ^ 1-4 
(O'Dell & Castaheda 1984). The observed range from galaxies at z ~ 1 in the DEEP2 survey is 
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similar, from 1.26 to 1.41, with a typical ratio of « 1.36 for the high equivalent width emitters (J. 
Newman, personal communication). If we fix the line ratio at 1.3 we find that simulated spectra 
yield a\(0 II) > 1.78 in only 6.7% of simulations. The corresponding probability for a/ is slightly 
higher at 12.9%. Thus, the best interpretation of this line is clearly Lyman-a, from both its 
asymmetry and its equivalent width. 

LALA J142441.20+353405.1: The second of these galaxies, LALA J142441.20+353405.1, was 
spectroscopically detected in both our Keck + DEIMOS data (covering 5160A to 9525A) and 
Gemini-N + GMOS data (covering 6177Ato 10484A). A strong emission line is seen at a wavelength 
9136A. This line is most likely [O III] A5007 at a redshift z = 0.824, based on a careful search for 
other emission lines. Unfortunately, the expected locations of the [O III] A4959 and H/3 lines both 
fall atop night sky OH emission features. The interpretation of this object depends critically on 
these features (with Lyman-a being the best interpretation if no secondary lines exist). 
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Fig. 3.— Spectrum of LALA J142441. 20+353405.1, obtained with the DEIMOS spectrograph on 
Keck II. The solid histogram shows the measured flux after smoothing with a 3 pixel median filter 
to suppress residual features from night sky lines. The dashed line shows the la flux uncertainties 
from photon counting errors. We have marked the wavelengths where the [O III] A4959 and H/3 
lines would fall if the primary line at 9136A is [O III] A5007. 
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Fig. 4.— 2 dimensional DEIMOS spectrum of LALA J142441. 20+353405.1. Labeled contours show 
the wavelength scale. Night sky lines are well subtracted here, but their locations are plainly visible 
as bands of enhanced photon noise. The presumed [O III] A4959 line is at 9048A, and is partly 
visible to the blue of the 9049A night sky line. The expected location of H/3 is at 8870A. 
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We therefore fitted Gaussian line profiles at the expected locations of both [O III] A4959 and 
H/3, using fitting weights derived from the observed profile of the 9136A line and from the observed 
noise as a function of wavelength in sky-subtracted spectra. This approach optimally exploits any 
part of the satellite line profile that is not obscured by night sky line residuals. We obtained a four 
sigma detection of the [O III] A4959 line and a line flux ratio of 0.33 ± 0.08 (/ (4959)// (5007)) from 
our 2-D DEIMOS data. A similar exercise using the 1-D GMOS data gave another 4cr detection 
at the expected [O III] A4959 wavelength and a line ratio of 0.35 ± 0.09. The consistency of these 
results with one another and with the theoretically expected line ratio of 0.33 leads us to believe 
this interpretation of the data. Further fitting yields also a marginal H/3 detection in the DEIMOS 
data (3a, line ratio 0.111 ± 0.04) but nothing significant in the GMOS data. Measurement of the 
[O II] AA3726,3729 lines with this method is more difficult, since the expected position of [O II] 
AA3726,3729 depends on the spectral trace (which we cannot measure because the continuum in 
LALA J142441. 20+353405.1 is undetected). In practice, we see a hint of the [O II] AA3726,3729 
doublet in both data sets, but it is more sig nificant in the GMOS data (2.4cr and 6Aa for 3726A 
and 3729A) than in the DEIMOS data (3.7cr and I.Oct). Further circumstantial support for the 
z = 0.824 interpretation comes from spectroscopic observations of a red galaxy sample in this field, 
which show a strong spike in the redshift distribution at z ~ 0.824 (Dey et al. 2004). 

Other foreground emission line interpretations are ruled out: An Ha line at 9136A should 
imply [O III] AA4959,5007 and H/3 near 6900A, which is not seen, and [O II] AA3726,3729 at 9136A 
would be easily resolved as a close doublet in our DEIMOS data. 

The flux of the 9136A line is large but subject to systematic uncertainties. Flux standards 
from our DEIMOS run imply a line flux of 3 < fn < 5.7, where fn = //(10~ 17 ergcm~ 2 s _1 ) is the 
scaled line flux. Our LALA survey narrowband imaging would naively imply a flux fn = 2.2, but 
the measured wavelength places the emission line on the wing of the narrowband filter, where the 
throughput is low. Correcting by a factor T max /T(9136A) implies a line flux 4 < fn < 6.6, where 
the lower value is based on the original filter specifications given to the vendor and the higher value 
is based on filter tracings made in May 2003 at NOAO by Heidi Schweiker and George Jacoby. We 
adopt fn = 4 as a reasonable best guess. 

The line is quite narrow, with an observed FWHM of 3.55A in our DEIMOS data (corre- 
sponding to 5.5 DEIMOS pixels, or a velocity width of ss 100km s _1 ). This width corresponds 
well to the typical seeing during our observations, which means that the physical line width in 
LALA J142441. 20+353405.1 could be < 100km s" 1 and also that the angular size of the line emit- 
ting region is no larger than 0.7arcsec. Applying our line asymmetry test here, we find C sym ~ 0.5. 
Because there is no demonstrable asymmetry, we did not pursue measurement of a\ and a/. A 
Gaussian fit to the line yields a x 2 /d.o.f. = 0.84 measured from 9132 to 9140A. 

We also see no evidence for continuum emission. The best equivalent width estimate for this 
source comes from our DEIMOS data, which yield a 2a lower bound of EW > 620A in the rest 
frame for z = 0.824. While large, this is within the observed range for [O III] A5007. We have 
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previously reported a z = 0.34 [O III] A5007 source with ~ 1000A rest equivalent width (~ 1400A 
observed), discovered by its narrowband excess in the LALA Bootes field (Rhoads et al 2000). 
Statistical studies of nearby HII galaxies show a range up to ~ 1700A (Raimann et al 2000). 

LALA J142441. 20+353405.1 and LALA J142442. 24+353400.2 lie 13" apart, corresponding to 
a projected separation of ~ 80kpc (physical). This proximity appears to be a coincidence given the 
differing line identifications. If the lines were the same transition, the velocity separation of the 
two would be 800km s _1 . 

Other candidates: Of the remaining two candidates, LALA J142544. 41+353327.7 was detected 
only in the 918nm filter and is probably a noise spike in that filter, based on its nondetection in 
Keck DEIMOS data. Assuming Gaussian statistics, we should expect ~ 1 noise peak above our 
5<7 detection threshold among the ~ 10 6 ' 5 independent resolution elements in the narrowband im- 
age, and more may be found if the noise properties of the image are not perfectly Gaussian. The 
second spectroscopic nondetection, LALA J142610. 55+354557. 6, may have been a time-variable 
continuum source rather than a true emission line object. It was targeted in a GMOS mask ob- 
served on 2003 July 02 and 03, with a total integration time exceeding that required to detect 
LALA J142441. 20+353405.1 and LALA J142442. 24+353400. 2, which have similar narrowband 
fluxes. It appears as a narrowband excess target when compared to images from 2001 and be- 
fore, but has a marginal R band detection in NDWFS images from UT 2002 April. This R band 
detection is offset from the 918nm source by about 0.8", and may or may not be the same source, 
since the faintness of the object renders precise narrowband astrometry difficult. Regardless of its 
precise nature, it appears not to be a redshift 6.5 Lyman-a galaxy. 

6. Discussion 

6.1. Physical Parameters of the Sources 

Adopting a cosmology with Q m = 0.27, $7 a = 0.73, and Hq = 71kms~ 1 Mpc~ 1 (cf. Spergel 
et al. 2003) gives a luminosity distance of 2.00 x 10 29 cm (65Gpc) for z s=s 6.53. We then ob- 
tain a line luminosity of (1.1 ± 0.2) x lO^ergs" 1 for LALA J142442. 24+353400.2. The other 
three Lyman-a galaxies so far spectroscopically confirmed at z ~ 6.5 have luminosities of 1.0 x 
10 43 ergs -i (g DF J132415.7+273058; K03), 0.56 x lO^ergs" 1 (SDF J132418.3+271455; K03), 
and 0.3 x lO^ergs" 1 (HCM 6A; H02). Thus, LALA J142442. 24+353400.2 is comparable to 
SDF J132415. 7+273058 and may be the most luminous Lyman-a line yet seen at this redshift. 
Comparison with Lyman-a luminosity function fits at z ~ 6.5 that we derive in a companion paper 
(Malhotra et al. 2004, in preparation) suggests that the observed line luminosity exceeds the char- 
acteristic Lyman-a luminosity L Lya,+ of emission line selected galaxies at this redshift by a factor of 
a few. While luminosities substantially above L* are often associated with strong gravitational lens- 
ing in high redshift source populations, we consider this unlikely for LALA J142442. 24+353400. 2, 
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since any massive lensing galaxy along the line of sight would have to have magnitude Iab ^ 25.8 
to avoid detection in the NDWFS I band image. 

We can convert the line luminosity to an estimated star formation rate using a conversion factor 
of lM /yr = 10 42 ergs -1 , which yields 11 M /yr. This conversion follows from a widely used set 
of assumptions: a Kennicutt (1983) initial mass function (IMF) and corresponding conversion 
between Ha luminosity and star formation rate, plus the Lyman-a to Ha ratio for dust-free case B 
recombination. These are probably not valid in detail for high redshift Lyman-a galaxies, but as we 
do not yet have a well justified model, it is convenient to at least use a standard one. Combined with 
our survey volume of 2.1 x 10 5 Mpc 3 (comoving), we obtain a lower bound on the star formation 
rate density (SFRD) of 5.2 x lO- 5 M yr- 1 Mpc -3 from LALA J142442. 24+353400.2 alone. Of 
course, given that LALA J142442. 24+353400. 2 is much brighter than L*, its contribution to the star 
formation rate density is only a small fraction of the total. For comparison, K03 probe further down 
the luminosity function and find a lower limit to the SFRD of 5 x lO~ 4 M yr _1 Mpc -3 , while the 
Lyman-a line found in the H02 survey would imply a la range from 3 x 10~ 3 to 6 x 10 -2 . All of these 
rates are subject to incompleteness corrections for sources falling below survey detection limits in 
either Lyman-a flux or equivalent width, and to additional systematic uncertainties associated with 
the stellar population model assumed and with the radiative transfer of the resonantly scattered 
Lyman-a line. In particular, the model used here and elsewhere to convert between Lyman-a 
luminosity and star formation rate cannot actually reproduce the observed Lyman-a equivalent 
width distribution at z ~ 4.5 (Malhotra & Rhoads 2002). Still, relative comparisons among Lyman- 
a samples remain valid, and show that sources like LALA J142442. 24+353400. 2 are rare, luminous, 
and constitute a fraction < 10~ 2 of the global Lyman-a production and star formation rate at 
z pa 6.5. 

For LALA J142441. 20+353405.1, we have flux measurements or limits on the [O III] A5007, 
[O II] A3727, and H/3 transitions. These allow us to place general constraints on the metallicity of 
this galaxy using R 23 = ([/(5007) + /(4959) + /(3727)]//(#/?)). We find log(i? 23 ) > log([/(5007) + 
/(4959)] //(#/?)) = 1.1 ± 0.12. The additional contribution to the oxygen flux from [O II] A3727 is 
either negligible (using the DEIMOS measurement, for which 32 = [/(5007) + /(4959)]//(3727) » 
20) or small (using the GMOS measurement, for which O32 ~ 2.5). From Kewley & Dopita (2002), 
the theoretical maximum for log(i?23) is 0.97, achieved only for a modestly sub-solar metallicity 
near \og(0/H) + 12 ~ 8.38 and a high ionization parameter q > 3 x 10 8 cm/s. (Here the ionization 
parameter is defined as the ratio of ionizing photon flux to number density of hydrogen; see Kewley 
and Dopita 2002.) Our observation is consistent with this theoretical maximum, but not with i?23 
substantially below this maximum, and we thus conclude that LALA J 142441. 20+353405.1 has 
12 + \og{0 /H) ~ 8.4. For comparison, estimates of the solar value for 12 + log(0/H) range from 
8.9 (McGaugh 1991) to 8.68 (Allende Prieto, Lambert, & Asplund 2001). 

This places LALA J142441. 20+353405.1 in the lowest quartile of metallicity for galaxies at 
0.47 < z < 0.92 from the Canada- France Redshift Survey sample (Lilly, Carollo, & Stockton 2003), 
but around the midpoint of the metallicity range for Lyman break galaxies at z « 3 (Pettini et 
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al 2001). The ratio O32 also shows a closer physical resemblance to LBGs than to continuum- 
bright galaxies at 0.47 < z < 0.92: Ratios O32 > 2 are observed in only one of 66 galaxies 
reported by Lilly et al (2003) but in four of five LBGs in the Pettini et al (2001) sample. However, 
the apparent discrepancy with other z ~ 0.8 galaxies is likely due to the larger characteristic 
luminosity of the Lilly et al sample, which has Iab < 22.5. Our I band limiting magnitude 
(Iab > 25.9) corresponds to Johnson Mb > —16.9. Using a metallicity-luminosity relation derived 
by Melbourne & Salzer (2002) for local emission line selected galaxies, we would expect 12 + 
log(0/H) = 8.1+0.27 where the error bar is the empirically determined scatter for the local sample. 
Thus, LALA J 142441. 20+353405.1 actually has a slightly high metallicity for its luminosity. This 
is broadly consistent with the finding of Lilly et al (2003), that most galaxies at 0.47 < z < 0.92 
have metallicity comparable to local galaxies of similar luminosity. 

Photoionization models show that equivalent widths of ~ 600A correspond to starburst ages 
of 3 to 4Myr for metallicity Z ~ Z Q /4 (Stasihska & Leitherer 1996), implying that the emission 
from LALA J142441. 20+353405.1 is dominated by a very young stellar population. Combining 
this age with the gas metallicity implies that this is not the first generation of stars formed in this 
high-redshift dwarf galaxy. The current star formation rate can be estimated from the emission 
line fluxes. The conversion between line flux and star formation rate is better established for H(3 
(using a ratio f(Ha)/f(H/3) = 2.8 and the Kennicutt (1983) conversion factor) than for [O III] 
A5007. The observed HP flux thereby implies a star formation rate is ~ O.4M /yr, though the 
combined uncertainty in the line flux and conversion factor is at least a factor of two. 



6.2. Implications for Reionization 

Because Lyman-a photons are resonantly scattered by neutral hydrogen, Lyman-a emitting 
galaxies may suffer considerable attenuation of their line flux when embedded in an intergalactic 
medium (IGM) with a substantial neutral fraction. Thus, a decrease in Lyman-a galaxy counts 
may be expected at redshifts substantially before the end of reionization (Rhoads & Malhotra 2001, 
Hu et al. 2002, Rhoads et al. 2003). Typical high redshift Lyman-a galaxies have small continuum 
fluxes, and upper limits can be placed on their expected Stromgren sphere radius r s by combining 
their observed fluxes and equivalent widths with stellar population models (Rhoads & Malhotra 
2001). A radius of r s ~ 0.5Mpc is typically inferred. For comparison, r s > 1.2Mpc is required to 
avoid a scattering optical depth r > 1 at the systemic velocity due to the neutral IGM outside the 
Stromgren sphere. 

Haiman (2002) has argued that the first reported z ~ 6.5 galaxy (Hu et al. 2002) could still 
be embedded in a neutral IGM, because (a) emission from the red wing of the emitted line will be 
less strongly scattered than line center, reducing the effective attenuation of the line, and (b) the 
observed Lyman-a equivalent width in this object is substantially below the theoretical maximum 
from stellar population models, so that the observed line could in fact be substantially attenuated. 
Santos (2003) considers additional physical effects, including gas motions and overdensity associated 
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with the formation of Lyman-a galaxies, which generally increase the Lyman-a attenuation over 
that in Haiman's model. For example, in the case of infall to a forming galaxy, scattering by residual 
neutrals within the Stromgren sphere can produce large optical depths (r ~ 100) on the blue side 
of the Lyman-a line, between the systemic velocity and the gas infall velocity. However, Santos 
also suggests that the velocity differences observed between Lyman-a lines and other emission lines 
in Lyman break galaxies could be intrinsic to the source, so that the entire Lyman-a line is emitted 
with an effective redshift of a few hundred km s _1 . In this case, the effective Lyman-a attenuation 
for z ~ 6.5 is a factor of ~ 3 ; and is not very sensitive to model parameters besides the ionization 
state of the IGM. 

LALA J142442. 24+353400. 2 constrains reionization better than most other z k, 6.5 galax- 
ies because of its relatively large equivalent width (> 70A rest frame, 2a, before corrections for 
IGM opacity; or > 40Awith corrections suitable for an ionized IGM; see section 5). The source 
SDF J132415. 7+273058 (K03) has a similar line luminosity, a continuum detection near the upper 
limit for LALA J142442. 24+353400. 2, and hence an equivalent width measurement quite similar 
to LALA J142442.24+353400.2's 2a limit. 

In a neutral universe, these sources would require an intrinsic equivalent width > 120A even 
if the Lyman-a line is redshifted as suggested by Santos, and > 400 A in the absence of such a 
redshift. The former is consistent with the observed Lyman-a equivalent width distribution at 
z ~ 4.5 (Malhotra & Rhoads 2002). The latter is high even for the z ~ 4.5 sample. Thus, the 
properties of these objects are most easily understood if the universe is mostly ionized at z ~ 6.5. 
More robust constraints on reionization from individual Lyman-a galaxies will be difficult, because 
the expected suppression of the Lyman-a line depends strongly on the velocity offset between 
emitted Lyman-a and the galaxy's systemic velocity, and measuring the latter would in general 
require spectroscopy in the thermally dominated mid-infrared (Santos 2003). 

Stronger constraints on reionization from Lyman-a galaxies will be possible using statistical 
samples. Suppression of Lyman-a by the IGM will have a strong effect on the Lyman-a luminosity 
function: Even a factor of ~ 3 reduction in Lyman-a flux will reduce the observed Lyman-a number 
counts dramatically. The accessible portion of the Lyman-a luminosity function is very steep (see 
Malhotra et al. 2004), and a factor of 3 in luminosity threshold corresponds to a factor of > 10 in 
expected Lyman-a source counts at fixed sensitivity. Moreover, the effect on the equivalent width 
distribution will be identical to that on the luminosity function, which would not be generically 
expected for other factors that could modify the Lyman-a luminosity function (such as evolution). 



7. Conclusions 

We have performed a narrowband search for Lyman-a emitting galaxies at z ~ 6.5 in the 
Bootes field of the Large Area Lyman Alpha survey, and have obtained spectra for all of our 
viable candidates. One source, LALA J142442. 24+353400. 2, is confirmed as Lyman-a at z = 
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6.535, based on an isolated, asymmetric emission line. The Lyman-a luminosity of this source is 
large, ps 1.1 x 10 43 ergs _1 . Objects this bright in the line are rare, occurring at a rate ~ 1 per 
2 x 10 5 comovingMpc 3 (this work and K03). It is therefore likely to correspond to a high peak in 
the density distribution at redshift z = 6.535. The equivalent width is also larger than most other 
z ps 6.5 galaxies, with a 2a lower bound of > 40A (rest frame, corrected for IGM absorption). This 
is difficult to reconcile with a neutral intergalactic medium unless the Lyman-a line is intrinsically 
strong and is emitted from its host galaxy with an intrinsic Doppler shift of several hundred km s _1 . 
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